function [next_pos] = ForageGWO(a, tar_pos, current_pos)

r1=rand();
r2=rand(); 
A1=2*a*r1-a; % Equation (3.3)
C1=2*r2; % Equation (3.4)
D_alpha = abs(C1*tar_pos(1,:)-current_pos(1,:));
next_pos(1,:) = tar_pos(1,:) - A1*D_alpha;

end